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Numerical search for a fundamental theory 

Vitaly VanchurirQ 

Arnold-Sommerfeld-Center for Theoretical Physics, Department fiir Physik, 
Ludwig-Maximilians-Universitat Miinchen, Theresienstr. 37, D-80333, Munich, Germany 

We propose a numerical test of fundamental physics based on the complexity measure of a general 
set of functions, which is directly related to the Kolmogorov (or algorithmic) complexity studied in 
mathematics and computer science. The analysis can be carried out for any scientific experiment 
and might lead to a better understanding of the underlying theory. From a cosmological perspective, 
the anthropic description of fundamental constants can be explicitly tested by our procedure. We 
perform a simple numerical search by analyzing two fundamental constants: the weak coupling 
constant and the Weinberg angle, and find that their values are rather atypical. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The search for a fundamental theory of nature has puz- 
zled many generations of scientists. Over time, mathe- 
matics describing physics became extremely complicated, 
but nevertheless very successful. Due to complexity, 
many branches of physics are evolving without much 
overlap simply because it is hard for a single human being 
to keep track of all of the new developments. Unfortu- 
nately, that might also mean that we could be missing 
something very fundamental that relates different fields, 
different theories, or even different constants of nature in 
a somewhat trivial way. 

In this article we intend to test if the constants of na- 
ture are truly independent. In Ref. 3] an overall twenty 
six fundamental constants are listed, with presumably 
another six required to fit the upcoming data. A natu- 
ral question to ask is whether all of these constants are 
really independent, or if some of them can be expressed 
in terms of the others. Ideally, one would hope to reduce 
the number of fundamental constants to just a few, the 
fewer, the better. 

If there are no fundamental reasons that fix some con- 
stant C to the measured value Cm, then we should be 
able to vary it and still obtain corresponding results. In 
other words, if we leave all of the physics unchanged, and 
add A to Cm, then the modified theory, with the constant 
given by C = Cm + A, should approach the original one 
in the limit A ^ 0. 

On the other hand, if the observed value of the con- 
stant Cm is described by an anthropic selection process 
P, then it must be obtained from some probability 
distribution P{C). If P{C) is smooth, then in a small 
enough neighborhood of Cm the function is nearly flat. 
This means that we are equally likely to observe an arbi- 
trary value of C within a small neighborhood of the mea- 
sured value Cm- Thus, the anthropic and non-anthropic 



views of fundamental physics seem to support the follow- 
ing hypothesis: A fundamental constant C must not lie 
in a preferred location on the interval {Cm — A, Cm + A) 
for small A. 

In the following sections, we will define exactly what 
we mean by preferred location. In contrast, if we find 
that some constant takes a highly atypical value withiir 
its neighborhood, then the constant in question may not 
be fundamental, or its anthropic function P{C) is highly 
irregular. 

The paper is organized as follows. In the second section 
we describe the basic idea by preforming a simple numer- 
ical search. In the third section we extend the analysis 
to an arbitrary pair of fundamental constants and apply 
the procedure to a simple model. In the forth section a 
general mathematical formalism is developed. In the last 
section we summarize the main results. 



II. MOTIVATIONS 

Measurements in any experiment could be described by 
a set of real numbers with some uncertainties^ . If some 
experiment is explained by a theory, then there are mea- 
surements that are related by expressions derived from 
the underlying theory. On the other hand, there are al- 
ways some other measurements that are completely inde- 
pendent of one another, based on the theory in question. 
We will refer to such measurements as fundamental con- 
stants of the theory. 



A. Numerical example 

To quickly motivate the analysis of this paper we con- 
sider two fundamental constants: the weak coupling con- 



^ All physical quantities are given in Planck units with c 
'Electronic address: ' vitaly@cosmos.phy.tufts.edu| G = 1. 
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stant g at mz and the Weinberg angle 9w- In dimen- 
sionless units the constants take the following values [1] : 



(a) Statistics with four eiementary functions 



,9 = 0.6520 ±0.0001 
0w = 0.48290 ±0.00005 



(1) 



We would like to perform a numerical test of the pre- 
cision intervals (0.6519,06521) and (0.48285,0.48295), 
to check weather these intervals are typical or weather 
they are "carefully chosen". (Intuitively, the in- 
tervals (3.1415926,3.1415928) or (0.9999999,1.0000001) 
seem to be very atypical, whereas the interval 
(4.2831741,4.2831743) is in some sense more random.) 

Let us choose only three constants tt, e, 1 and five ele- 
mentary functions: 



fi{x,y) = x + y, 
h{x,y) ^ x-y, 
f3{x,y) = xy, 
fi{x,y) = x/y, 
fb{x) = \/x. 



(2) 



For the time being we define the complexity measure M 
of expressions as the total number of elementary func- 
tions in the expression. For example. 



M\ 



(1 + 1) + - 



■e- ^ (e-e) - 1 ' 

M[/4(/i(/i(l, 1), tt), /2(/3(e, e), 1))] = 5. 



(3) 



The numerical complexity of an interval is then defined as 
the smallest complexity of all expressions that generate 
a number in the interval. For example, the complexity 
of the interval (0.9999999, 1.0000001) is 0, and the com- 
plexity of (1.9999999,2.0000001) is Af[/i(l, 1)] = 1. 

We have developed a computer program which can 
search through all possible expressions with complex- 
ity up to 6. The numerical complexities of 300 neigh- 
boring intervals around the measured value of 9w are 
plotted on Fig. [T] . All five elementary functions of 
Eq. ^ are used to generate the plot on Fig. [T]d, and 
only the first four functions are considered for the 
plot on Fig. [T^. In both cases, it appears that the 
numerical complexity of the interval (0.48285, 0.48295) 
is way too small, when compared to other inter- 
vals: ... , (0.48265,0.48275), (0.48275,0.48285), 
(0.48295,0.48305), (0.48305,0.48315), ... . In fact, the 
Weinberg angle is given by a simple expression: 



1 



(4) 



with complexity measure equal to 3. 

The probability of observing such a small value of the 
complexity in the neighborhood of 9w is about 1.7% 
for the analysis with four elementary functions. (It 
could be seen directly from the Fig. [T^: only five out 
of 300 neighboring intervals have the complexity equal 




FIG. 1: Numerical complexities of 300 intervals around the 
measured value of the Weinberg angle. 



(a) Statistics with four eiementary functions 




0.63 0.64 0.65 

(b) Statistics with five eiementary functions 




O.'esV 0.66 
^ g=0.6520 



FIG. 2: Numerical complexities of 300 intervals around the 
measured value of the weak coupling constant ai mz- 



to 3.) The probability of observing the complexity of 
(0.48285, 0.48295) in the five functions analysis of Fig. 
[Hd is a bit higher and equal to 4.7%. Of course, that 
could have happened by chance, and thus one needs a 
better measurement of the Weinberg angle in order to 
prove, or disprove the above relation. 

The numerical complexity of the precision interval of 
g seems to be very typical in its neighborhood from the 
analysis with four functions (Fig. [2]a), and very atyp- 
ical (the probability is 2%) from the analysis with five 
elementary functions (Fig.^s). The following simple ex- 
pression was revealed: 
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(5) 



Once again, in order to prove or disprove the relation a 
much more precise measurement of g must be available. 



B. Discussion 

The reader must be wondering if we really propose that 
g and Qw are given by our simple expressions. The short 



answer is no. We have chosen the weak coupHng constant 
and the Weinberg angle only because out of all constants 
listed in Ref. they are measured with the highest 
precision. Nevertheless, the uncertainty was too large to 
draw a definite conclusion. Presumably, with better tech- 
nology we can hope to get much better measurements of 
the constants. It should also be noted that the precision 
could be dramatically improved by simply repeating the 
experiment. This is true even if the uncertainty princi- 
ple of quantum mechanics does not immediately allow to 
measure some constants precisely. If the precision of a 
single measurement is S, the combined uncertainty of N 
measurements is given by S/^/N, and thus can be arbi- 
trary small. 

Unfortunately, we were not able to study the proper- 
ties of the fine structure constant, which is already mea- 
sured with very high precision: a = 0.007297352568(24) 
. Not even a single expression with M < 7 was found 
to generate a number inside of the precision interval of 
a. To calculate the relation function of such a tiny in- 
terval, in general, the functions f £ G with much larger 
complexities must be checked, which is currently out of 
reach. Perhaps more advanced algorithms, such as PSLQ 
0, should be used to speed up the munerical search. 

Similar parameter fitting was first performed in Ref. 
0, where the expression a = 1/136 was claimed. The 
major difference in here, is that the statistical significance 
of such claims is actually verified. In fact, the "discovery" 
of Ref. could have been easily checked for a given 
set of elementary functions. Fortunately, the expression 
a = 1/136 (and a = 1/137) is already ruled out and we 
do not have to perform the analysis. However, there are 
other numerical coincidences that have been discussed in 
the literature. For example, in Ref. @], the Weinberg 
angle is given by the expression sin{9w) — \ ~ J2 '^tiich 
has a relatively large complexity. This should not be 
regarded as atypical, when the expressions with much 
smaller complexities could also describe the constant in 
question. 



III. RELATION FUNCTIONS 



In the previous section we have analyzed the numerical 
structure of some of the fundamental constants through 
the complexity of their precision intervals. The proposed 
procedure enabled us to search for the numerical expres- 
sions that could explain the observed values of these con- 
stants. In addition, we were able to quantify the statis- 
tical significance of the obtained simple expressions. In 
this section we will generalize the approach to study the 
numerical structure of an arbitrary pair of fundamental 
constants. 
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FIG. 3: Expectation value of the two-interval relation func- 
tion is calculated by averaging over all pairs of 5 intervals 
inside A-neighborhoods of Ci and C2. 

A. Two intervals example 

For simplicity, let us consider two fundamental con- 
stants given by positive real numbers Ci and C2 with 
precisions 5i/2 — 82/2 = 5/2. We assume that 5 <C 
|Ci| < IC2I or, in other words, Ci and C2 are measured 
very precisely. In what follows, we will discuss a sim- 
ple procedure which could give a definite answer to the 
question of whether constants Ci and C2 are really inde- 
pendent. 

We define a simple set of real valued functions and a 
measure defined on them. Let T be the set of functions 
of the form f{x) = mx and their inverses f^^{x) = x/m, 
where m € N. The measure M is a map from J- to N, 
defined as M{f) = m if f{x) = mx or f{x) = x/m. 

We can now define a two-interval (/i and I2) relation^ 
function i?(/i,/2) with respect to the set of functions J- 
and measure M: 

i?(/i,/2) = min{M(/i) + M(/2)| 

3zeM:/i(z)e/iA/2(z)e/2}. ^ ' 

(In other words, we are trying to approximate the ratio 
C1/C2 with two integers m,i/m2 ~ C1/C2, such that 
TOi -|- 7712 is the smallest.) 

In what follows, the two-interval relation function of 
the uncertainty intervals of Ci and C2 is denoted as: 

^0 ^ nm - + iC2 -\.C2 + -^\. (7) 

To calculate an expectation value of i?o, we have to av- 
erage the relation function over all pairs of 8 intervals 
inside larger neighborhoods (A ^ S) of C\ and C2 (See 
Fig. [3] ). In general, a numerical search must be per- 
formed in order to estimate the expectation value, but 
for our simple example we can find it analytically. 

The aim is to find two integers m\ and m2 correspond- 
ing to functions j\[x) — m\x and /2(a;) = 77122;, such that 
for some a; G M: 

|Ci-a;TOi|<f 

\C2-xm2\ < f ^ ^ 



^ Unlike the correlation function, which is defined on a set of 
points, the relation function is defined on a set of intervals. 
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In order to minimize M(/i) + M{f2) = mi + m2, it is 
enough to find the smallest m2, for which both inequal- 
ities are satisfied (for some mi). Since C2 > Ci, then 
TO2 > 'Til, and the two inequalities can be reduced to a 
single one; 



m2 



mi\ < S. 



(9) 



For any fixed m2, the probability that the inequality 
in Eq. ([9]) is satisfied (for some nii) can be written as 



P{m2) = 



11126 



(10) 



Thus, one would expect the inequality to hold for at least 
one value of 7712 if 



fe=i 



Co 



or 




"^2 ;C A/ -r 



(11) 



(12) 



This immediately implies that the expectation value of 
the relation function is approximately given by: 



E 



1712 + mi 



m2 + r7i2 



Ci Ci + C2 
C2 \/5C'). 



(13) 



For a particular choice of Ci , C2 and 5 there are three 
possible outcomes of the procedure: i?o ^ E, Rq <^ E, 
E ^ Rq, which we will discuss separately. 



B. Interpretations 

Naturally, one would expect Rq ^ E, if constants Ci 
and C'2 are not simply related with respect to the set of 
functions J-' and measure M. However, even then it is 
still possible that Ci and C2 are related by some other 
set of functions with respect to some other measure. The 
procedure for generating the most general set of func- 
tions is rather complicated and will be discussed in the 
following section. 

By far, the most important possible outcome of our 
procedure is the discovery that Rq <C E. In this case, it 
is clear that we have reached a contradiction and some- 
thing is completely wrong with our understanding of the 
theory. The numbers that we thought are not related are 
in fact related by some very simple relation. After this 
discovery, one can go ahead and study the functional rela- 
tion between Ci and C2 , in order to reformulate the the- 
ory using this hint. (For example, Mendeleev's periodic 
table of elements could have been discovered by apply- 
ing the discussed procedure with exactly the complexity 
measure on the set of functions that we have introduced 
in this section The atomic masses of elements are 



simply related to each other with respect to and mea- 
sure M, and thus Rq ^ E.) 

Another interesting result of the procedure is some- 
what more controversial. E Rq implies that Ci and 
C2 are atypically unrelated (or anti-related). It means 
that they are "placed" in specific locations on the real 
line such that one cannot simply relate them by a given 
set of functions. This might be a direct evidence that 
our set of functions and/or the measure procedure is not 
suited to study the particular pair of constants. Alterna- 
tively, it might also mean that one has to introduce yet 
another (more) fundamental constant, that would make 
measured values of Ci and C2 more typical in their neigh- 
borhoods. 



IV. GENERAL ANALYSIS 

Let us denote the general set of functions by G- Ev- 
ery element of the set / is a function constructed from 
a finite combination of its fundamental members: num- 
bers (a finite number of numbers) , functions (elementary 
operations, exponential function and special functions), 
loops (unevaluated integrals, sums and products), recur- 
sive relation and inverses of all of the above. Note that 
the set of fundamental elements is finite, when the set of 
general functions is infinite, but countable. Each func- 
tion / of the set Q can be represented as a binary tree 
with functions associated with nodes (one or few) of the 
tree. The complexity M of such tree can, for example, 
be defined as the sum of complexities of all nodes. 

With respect to Q and some complexity measure M we 
can construct a two-interval relation function R: 

Riluh) = min{M(/)|3zi,Z2 G R : 

(5i 69 (14) 

- Cil < ^ A |Z2 - C2I < y A /(zi, Z2) = 0}. 

Similarly, we can generate an n-interval relation function: 

R{Ii, ...,/„) = min{A/(/)| 3zi...z„ e R : 

6, S (15) 

\zi -Ci\< ^...|z„ - C„| < ^ A /(zi...z„) = 0}, ' 

where n S N. Altogether, there are 2" — 1 independent 
relation functions with respect to the precision intervals. 
If all n constants are really independent, then all of the 
relation functions i?2, ■••^2"-i of uncertainty inter- 
vals (xi — Y , -I- Y ) should be of the same order as their 
respective expectation values Ei, E2, ...i?2'>-i. These val- 
ues are calculated by averaging over all intervals: 

li{di) = {x^ + [di - ^)(5j, a;^ + [d, + ]^)5i) ^^^^ 

C {xi - Ai,x, + Ai), 

where di G Z. For example, the expectation value of the 
relation function for all n constants is given by: 



i?2"-l — 



-y- 



i?(/i(di),...,/„(d„)), (17) 
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where N is the normahzation constant. 

In order to compare the values of relation functions, 
intervals of the same size must be used. The relation 
function is larger for a smaller interval, as we have seen 
in the previous sections (Eq. (fT3|) ). If we take the size of 
an interval to zero, the relation function either converges 
to a finite number, or diverges to infinity. In the former 
case the respective constants are proved to be related, 
whereas in the latter case the constants are proved to 
be independent. Of course physically it is not possible 
to perform a measurement with infinite precision, and 
therefore we can only prove the dependence or indepen- 
dence to a certain degree. 

To study the structure of a set of constants wc might 
need to know more than just expectation values of the 
relation functions. It is straightforward to define a prob- 
ability distribution Pk of relation functions in A.^ neigh- 
borhoods of the measured values of constants C^'s. For 
a given probability distribution Pk{R), we can calculate 
the probability of observing a measured value of Rk or 
smaller: 



Pk{R)dR. 



(18) 



If S'fc <C 1 or 1 — S'fc <C 1, then the constants in question 
are atypically related or anti-related respectively. How- 
ever, we have to be very careful, because for a large num- 
ber of relation functions a very atypical situation may 
occur according to the central limit theorem. 

In order to draw a definite conclusion, we have to aver- 
age separately the values of Sk in the ranges S < ^ and 



n 



A = 



Er=T'(i-^fc)e(5fe-i) 



where Q{x) is a step function: 

ew4 ? 



if X < 
if X > 0. 



(19) 



(20) 



(21) 



Strictly speaking, TZ and A are only defined if the denom- 
inators are non-zero. This is the case for a large number 
of fundamental constants, and that is exactly when these 
expressions are particularly useful. We call TZ and A the 
relation and anti-relation indexes of a given set of con- 
stants. 

Some interpretations of different outcomes of our pro- 
cedure on a set of fundamental constants were given 
above. We can summarize the possibilities in the lan- 
guage of relation and anti-relation indexes. If the con- 
stants are: 

- independent, then TZ ^ A 1., 



- related, then TZ -^1 and ^ ~ 1, 

- anti-related, then TZ ^ 1 and A<S^1, 

- related and anti-related, then 7?. ^ 1 and ^ ^ 1. 

In the above analysis we have assumed that all relation 
functions are uncorrelated. If correlations between some 
relation functions are not negligible, the definitions for 
relation and anti-relation indexes have to be adjusted. 



V. CONCLUSION 

In this article we have developed a procedure of testing 
whether the fundamental constants are really fundamen- 
tal. The method relies on the hypothesis that the mea- 
sured value of a fundamental constant must not lie in 
a preferred location within its neighborhood. If the hy- 
pothesis is proved to be wrong, then the anthropic func- 
tion (if exists) must be sharply peaked around the mea- 
sured value of the constant. Similarly, a non- anthropic 
point of view implies that some constants of nature are 
very carefully chosen. Both approaches suggest that one 
might be able to derive the values of fundamental con- 
stants from first principles. 

As we have discussed above, our analysis relies on the 
definition of the complexity which is not unique. Never- 
theless, there are statements that one can make regard- 
less of the exact definition, given that all of the finite 
expressions have a finite complexity. In the limit when 
the size of the interval goes to zero, its complexity grows 
to infinity. This is true for all of the intervals, unless the 
interval must contain (due to some fundamental reason) a 
constant which is given by some fundamental expression. 
If the expression exists, then for very precise measure- 
ments of the fundamental constants it is guaranteed to 
be discovered. (In this case, an expectation value of the 
relation function can grow without bounds, whereas the 
relation function of the intervals in question is bounded 
from above.) 

Although wc have shown that in principle the exact 
definition of the complexity is not crucial, in reality the 
complexity must be defined such that a useful result is 
obtained as soon as possible. Thus, it is important to 
carefully choose the complexity measure, which would 
unable us to tackle the problem efficiently. Fortunately, 
similar ideas had been already proposed in mathematics 
and computer science 0, [l^, [lH. The algorithmic (or 
Kolmogorov) complexity is defined as a measure of the 
computational resource (e.g. the length of the shortest 
program) needed to specify an object. The relation of the 
complexity of intervals to the Kolmogorov complexity is 
an important question to be discussed in the upcoming 
publication. 

The last point that we would like to address is not com- 
pletely settled. The question is whether in practice it is 
possible to use our procedure to check widely accepted 
theories and/or to discover new physics. The discussed 
procedure relies very strongly on two parameters: preci- 
sion of measurements and computational resources. One 
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has to perform very precise measurements of fundamen- 
tal constants and then run the computer program until a 
meaningful result is produced. Both of these parameters 
depend very strongly on the state of available technolog- 
ical resources. However, even if we will discover that it is 
not possible to perform our analysis at present, this will 
undoubtedly change in the future, perhaps with help of 
quantum computers. 
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